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Abstract. This paper is devoted to investigating non-equiUbrium phase 
transitions to an absorbing state, which are generically encountered in reaction- 
diffusion processes. It is a review, based on [1, 2, 3], of recent progress in this field 
that has been allowed by a non-perturbative renormalisation group approach. We 
mainly focus on branching and annihilating random walks and show that their 
critical properties strongly rely on non-perturbative features and that hence the 
use of a non-perturbative method turns out to be crucial to get a correct picture 
of the physics of these models. 



1. Introduction 

Our understanding of equilibrium critical phenomena has largely benefited from the 
success of perturbative renormalisation group (RG) techniques and from conformal 
symmetry properties in two dimensions. Unfortunately, this success story has not 
spread out to non-equilibrium systems, though the latter are far more common in 
nature. Indeed, the study of phase transitions between non-equilibrium steady states 
has taught us that critical phenomena or generic scale invariancc turn out to be much 
richer far from thermal equilibrium — where detailed balance relations are violated 
— than in equilibrium statics or even near-equilibrium dynamics. 

However, despite considerable efforts, the very ingredients settling the universality 
classes out of equilibrium remain on a fragile footing. The problem originates mostly in 
the absence of an effective free energy functional that would allow to straightforwardly 
classify the universal behaviours in terms of symmetries and interactions. And even 
when such a functional exists — which occurs for reaction-diffusion processes or for 
Langevin type dynamics — analytical progress turns out to be difficult. On the one 
hand, models with Langevin dynamics cannot be conformal invariant. On the other 
hand, the efficiency of standard RG approaches is hindered by the fact that first, 
critical dimensions, when they are identified, happen to lie far from the physically 
interesting ones and furthermore, no calculations are available above two-loop order, 
which prevent from resorting to powerful re-summation techniques. In this context, 
any theoretical tool that can overcome the previous pitfalls is valuable and the non- 
perturbative renormalisation group (NPRG) stands as a promising candidate [4]. 
Indeed, this method has lead to great successes for systems at equilibrium during 
the last decade, appearing as a well-adapted tool to tackle strong coupling problems 
[4, 5] . Moreover, the NPRG formalism has been recently extended to non-equilibrium 
systems [1] and has given rise to important progress in the field of reaction-diffusion 
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processes, unveiling genuinely non-perturbative effects [1, 2, 3], which we review in 
this paper. 

This paper is organised as follows. In section 2, we give an overview of reaction- 
diffusion processes and recall their field theoretical formulation. In section 3, we briefly 
introduce the NPRG formalism generalised to non-equilibrium systems and wc derive 
the flow equations related to reaction-diffusion processes, before focusing on specific 
models. Some results on the universal properties of directed percolation (DP) are 
first reviewed in section 4 and the remainder of the paper (section 5) is devoted to 
the study of branching and annihilating random walks (BARW), which split into two 
categories, namely the 'odd' and 'even' BARW. 

2. Reaction-diffusion processes 

Reaction-diffusion processes constitute simple models that allow to gain some insights 
in non-equilibrium critical phenomena [6]. They are defined by a set of microscopic 
rules that govern the dynamics of a species of particles A. These particles randomly 
diffuse at a rate D and undergo some reactions, typically birth A (m -|- 1)A and 

death kA processes. From these competing interactions emerges at long time 
a stationary state which nature depends on the reaction rates Afe and am- Either all 
particles eventually vanish, leaving an empty state where all stochastic fluctuations 
cease and which is therefore called 'absorbing'. Or the density eventually saturates to a 
finite average value yielding an 'active' state where the dynamics constantly generates 
density fluctuations. The active and absorbing states are separated by a continuous 
phase transition [6]. 

A large range of the absorbing transitions encountered in reaction-diffusion 
processes fall into the DP universality class [7], which stands as the most prominent 
one. It has lead to a famous conjecture in the early eighties by Janssen and 
Grassberger [8] stating that a continuous transition to an absorbing state driven by 
a one-component order parameter will generically fall into the DP universality class 
(provided there is no additional symmetry or quenched disorder). The DP model, 
which corresponds to the rules 

A^2A, 2A^0, A^0, (1) 

therefore plays a paradigmatic role, as the counterpart of the Ising model for equilib- 
rium statistical physics. 

Reaction-diffusion processes naturally lend themselves to Monte-Carlo simula- 
tions, which have indeed largely contributed to our understanding of these processes 
(see [6, 9] for reviews). On the other hand, the simplest analytical approach is to de- 
vice a rate equation for the time-dependent average density n{t), assuming the various 
reactions to occur proportionally to the concentration of reactants. For instance, for 
the DP processes (1), this yields: 

dtn{t) ^ {a - ti)n{t) -2Xn(tf. (2) 

Eq. (2) entails a mean-field type of approximation since the density correlations are 
neglected — the joint probability of finding two particles at the same position has 
been factored. Eq. (2) is readily solved. It exhibits two stationary solutions Ua = 
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and rig = (a — /u)/(2 A), which stabihty depends on the sign of the 'mass' A = {cr — 
If A < 0, the density decreases to the absorbing solution n^, whereas it satmates 
to the active solution for A > 0. The explicit solution for the time-dependent 
density n{t) = nons/[no + {us — no) exp(— At)] shows that both asymptotic states are 
reached exponentially fast in time. The relaxation time A"-'^ diverges when a — fj, 
causing an algebraic decay of the density which corresponds to the critical state. 
The DP absorbing phase transition is illustrated in Fig. 1 [6], which represents the 
time evolution of particles on a one-dimensional lattice for increasing values of A (see 
caption) . 



Figure 1. Time evolution of particles on a one-dimensional lattice, starting from 
either a uniformly distributed initial configuration (top row) or a seed particle 
(bottom row), for the DP model (endowed with the dynamical rules (1)) from 
[6]. In this model, the stationary states (respectively absorbing on the left-hand 
side and active on the right-hand side) are reached exponentially fast in time, 
and arc separated by a critical state (middle column) where both time and space 
correlation lengths diverge, resulting in an algebraic decay of the density. 




Absorbing phase transitions can be characterised by a set of critical exponents, 
typically: 

ris-^ip-Pcf (3) 
U-^\P-Pcr^ (4) 

^11 ~ \P-Pcr^ (5) 

where z = i'\\/i'± embodies the dynamical exponent, which represents the anomalous 
scaling between space and time. 

From rate equations such as (2), one can simply work out the — d-independent 
— values of the critical exponents at mean-field level (for instance for DP /3 = 1, 
uj_ = 1/2 and 2 = 2 [6]). 

Field theory 

The mean-field picture holds as long as the density remains homogeneous — short- 
range correlated — enough in the system so that the role of space and time 
fluctuations is indeed negligible, which is generally justified in high dimensions. 
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However, the reaction-diffusion processes under scrutiny are so-called 'diffusion- 
limited', which means that the diffusion is not efficient enough to compensate the effect 
of reactions that locally alter the density distribution and hence invalidate the mean- 
field approximation. A finer analysis therefore requires a systematic incorporation 
of spacio-temporal fluctuations, which can be achieved through the construction of a 
field theory. 

A field theory to describe reaction-diffusion processes can be derived following 

a well-known formalism due to Doi and Pcliti [10]. The starting point is the 
master equation, which describes the time evolution of the probability distribution 
P{{ni},t) of the occupation numbers Ui of the sites i of a lattice — assuming no 
occupation restriction. The idea is to write the change in the occupation numbers 
on each site by means of creation and annihilation operators. Upon introducing Fock 
states \{ni}) to represent the configurations of the lattice, and a formal state vector 
I'^'l^)) = S the master equation can be recast into a Schrodinger-like 

equation for the state vector dt\^{t)) = —H\^{t)). Then, resorting to coherent- 
state path integrals as in quantum mechanics, one can construct a functional integral 
(p] = J 'D(f)'Di(f)Cx:p{—S\(j), 0]), which captures exactly (up to the continuum space 
limit) the microscopic stochastic fluctuations [10, 11]. The time-dependent statistical 
averages of observables — which are necessarily mere functions of the occupation 
numbers — can then be computed from Z. For the general processes A {m-\- 1)A 

and kA 0, this procedure yields [12]: 



The diffusion is encoded in the kinetic part and stands as the Gaussian (quadratic) 

theory which corresponds to the Brownian motion. All the reactions give rise to 
potential interaction terms. The action (6) can then root field theoretical treatments, 
and in particular NPRG methods. This action indeed underlies in the following the 
investigation of various reaction-diffusion processes, namely DP and BARW models. 
Before presenting these analyses (in sections 4 and 5), we first give a brief overview of 
the NPRG techniques. 

3. Non-perturbative renormalisation group out of equilibrium 

We do not intend to detail here the implementation of the NPRG, but rather set out 
its principle (and refer to [5, 4] for reviews). This formalism relies on Wilson's RG 
idea [13], which consists in building a sequence of scale-dependent effective models, 
that interpolate smoothly between the short-distance physics at the (microscopic) 
scale A; = A and the long-distance physics at the scale = 0, through progressively 
averaging over fluctuations. Rather than expressing — as in the original Wilsonian 
formulation — the flow of effective Hamiltonians for the slow modes, one can work out 
the flow of effective 'free energies' for the rapid ones, following [14, 5]. thus only 
includes fluctuation modes with momenta \q\ > k. At the scale fc = A, no fluctuation 
is yet taken into account and Fa coincides with the microscopic action S, while at 
fc = 0, all fluctuations are integrated out and Fq is the analogue of the Gibbs free 
energy F at thermal equilibrium, in that it encompasses the long distance properties 
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of the system. To construct T^, one needs to freeze the slow modes. This is achieved 
by adding a scale-dependent mass-like term to the original action [5, 1, 15]: 

asSA\= I 0(x,^)i^fe(v^5t)</>(x,^), (7) 

where the cutoff function Rk behaves as Rkiq^jUj) ~ fc^ (in Fourier space) for small 
momenta \q\ < k — so that the slow modes are decoupled — and Rk vanishes for large 
momenta \q\ > k — so that the rapid ones remain unaltered. The "partition functions" 
^k\j,j] = / 'D^Vi^ exp(— <S — A<Sfe -|- J j <P + J j ^) become therefore fc-dependent. Ffe 
is obtained through the Legendre transform of log Zk [j, j] : 

Tfe log Zk [j, j] = Jj^lj + J U - A5fe [V', ^] (8) 

and is a functional of the conjugate fields tp = 5\ogZk/5j and ip = S\ogZk/Sj . Note 
that the last term in Eq. (8) ensures that Fk has the proper limit at /c = A: Tk=A ~ S 
[5, 15]. An exact functional differential equation governs the RG flow of Ffe under an 
inflnitesimal change of the scale s = log(fc/A) [5, 1]: 

d,Tk = ^Tr / dA (f + Rk) , (9) 

where Rk is the symmetric, off-diagonal, 2x2 matrix of element Rk (in Fourier space) 

and FJj, [xp, xp] the 2x2 matrix of second derivatives of F^ with respect to ip and ip. 
Obviously, Eq. (9) cannot be solved exactly and one usually devices an ansatz for F^. 
However, as the approximations used do not rely on the smallness of a parameter, 
the approach remains non-perturbative in essence. In particiilar. it is not conflned to 
weak-coupling regimes or to the vicinity of critical dimensions and is therefore suitable 
to overcome the usual perturbative RG schemes. 

Since the critical physics corresponds to the long distance {q 0) and long time 
{oj 0) limit, a sensible truncation consists in expanding Tk in powers of gradients 
[14] and time derivatives. Retaining only the leading order in derivatives, the most 
general ansatz for F^ related to the field theory (6) writes [1]: 

Tk{i^,i>) = J d'^xdt{i>{x,t)[Dk[iP,i>]dt-Zk[iPj]V^]i^ + Uk[i'J]. (10) 

The effective potential Uk encompasses all possible reactions. The kinetic 
renormalisation functions Dk and Zk of the diffusive part account for the anomalous 
scalings of the fields and of time. Indeed, the anomalous dimension 77 of the fields and 
the dynamic exponent z are defined such that at criticality, tptp ~ k'^^^ and iv ^ k^ 
respectively [16]. It follows that r]k = —dglnDk and Zk = 2 + r]k + dslnZk, the critical 
exponents r/ and z corresponding to the (fc-independent) fixed point values of rik and 

Zk- 

The flow equations for the potential Uk and the renormalisation functions Zk 
and Dk have been established for generic reaction-diffusion processes in [1]. Different 
levels of approximation can be implemented. The simplest one, the so-called local 
potential approximation (LPA), consists in neglecting the kinetic renormalisations. 
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i.e. in setting Zf. = Df. = 1, upon which z = 2 and r] = 0. The LPA generally enables 
one to get a reliable qualitative picture of the physics as well as a fairly accurate 
determination of the static exponent v. This approximation can be refined by including 
(field- independent) kinetic renormalisation coefficients Dk and Zk, which allows for 
non-trivial estimates of ry and z. This approximation is referred to as leading order 
(LO). Finally, to get more accurate values of the exponents requires to incorporate 
the field dependence of the kinetic functions Zk[tp,x(;] and Dk[ip,i'] — referred to as 
next to leading order (NLO) — which becomes much more tedious numerically. 

In order to study a specific model, the generic flow equations dsUk, dgZk and 
dsDk must be solved for the functional Tk respecting the symmetries of the problem, 
starting at scale A from the corresponding microscopic action (6). We focus in the 
next section on the DP universality class. 

4. Directed Percolation 

The action corresponding to the DP processes (1) can be deduced from (6) and writes: 

Jx,t 

It can be conveniently rewritten expanding the response field around its stationary 
solution (j){x, t) = 1 + 4>{x, t) and rescaling the fields according to ^ — > ^JlXja 4> and 

— > -\/a72A (j), which leads to: 

Sdp[<P,^]^ jd'^xdt\^4>[dt-Dy^-{a-ii)](j3+V2^[<P(P^-4>^(i3]+X{(j3$fy (12) 

This action turns out to be equivalent to the so-called Reggeon field theory [17], 
which has been studied in particle physics since the early seventies [18]. The critical 
exponents have been computed in that context to two-loop order. However, the upper 
critical dimension of this field theory is dc = 4|, which is distant from the physical 
dimensions d = 1 and d = 2. Moreover, in spite of its simplicity, the DP model has no 
exact solution in d = 1 — at variance with the Ising model. Thus, the best estimates 
of the critical exponents of DP rely on numerical calculations and are given in table 
1. 

Providing analytical estimates of the DP exponents thus requires a method that 

is not confined to the vicinity of a critical dimension. This motivates the use of a 
non-perturbative approach — the NPRG — to fulfil this task, since the latter allows 
to span arbitrary dimensions. The generic flow equations introduced in section 3 can 
be used to compute the critical exponents of the DP universality class, provided the 
symmetries of the DP model are specified. The bare action (12) turns out to be 
invariant under the change: 

f <P{x,t) -(p{x,-t) 
\ 4>{x,t) -(l){x,-t), 

which is called the 'rapidity' symmetry. The effective potential Uk{'ipi''P) and kinetic 
renormalisation functions Dk{tp,tl>) and Zk{il>,'ip) — denoted Xk in the following — 

1 Note that in this paper, d denotes the dimension of space — i. e. does not include the time direction. 



(13) 
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Table 1. Critical exponents of the DP universality class, from NPRG calcula- 
tions within different levels of approximation LPA, LO in d = 1,2,3 [1] and 
NLO in d = 2,3 [15]. The last column gathers the best numerical estimates, 
ensuing from Monte Carlo (MC) simulations [19, 20] and series expansions [21]. 



d 




LPA [1] 


LO[l] 


NLO [15] 


MC+serics [19, 20, 21] 






0.584 


0.548 


0.59 


0.581(5) 


3 




0.872 


0.782 


0.83 


0.81(1) 




z 


2 


1.909 


1.90 


1.90(1) 




V 


0.730 


0.623 


0.73 


0.734(4) 


2 


3 


0.730 


0.597 


0.59 


0.584(4) 




z 


2 


1.884 


1.70 


1.76(3) 




V 


1.056 


0.888 




1.096854(4) 


1 


p 


0.528 


0.505 




0.276486(8) 




z 


2 


1.899 




1.580745(10) 



must be invariant under the rapidity transformation (13). This in turn imposes that 
the generic term of the Taylor expansion of the Xk^s be of the form aai3{4>"ip^ + 
{—l)'^'^^'tp"ip^), which only involves the two invariant combinations x = iptp and 
y = tp — tp. Thus, parameterising the X^s in terms of the invariants x and y ensures 
that they satisfy the rapidity symmetry constraints. 

The flow equations for Uk , Dk and Zk then have to be solved numerically. One can 
integrate the flows with the scale s starting from some initial bare rates [1]. For a flne 
tuned initial mass A, the (dimensionless) effective potential flows to a fixed function, 
which corresponds to criticality. The exponents are calculated in the vicinity of this 
fixed solution. This procedure has been performed at different levels of approximation, 
namely the LPA, the LO and the NLO [1, 15]. The results are displayed in table 1 
and show that the estimates converge fairly rapidly as the approximation is enriched 
to values in good agreement with the best numerical estimates. 

5. Branching and annihilating random walks 

Initially introduced by Bramson and Gray [22], models of BARW can be seen as 
reaction-diffusion processes endowed with the generic rules: 

A-^ {m+l)A, kA^0, (14) 

from which is excluded the spontaneous decay A ^ {i.e. k > 1). The particles 
can therefore only disappear through a fc-body annihilation. This restriction has 
drastic consequences since at mean-field level, it forces the system to reach an active 
state with mean density = [{mam)/{k Xk)Y^'^'^"^^ for any non-zero branching rate 
dm- A trivial transition occurs at am = 0, where the model coincides with the pure 
annihilation model. The latter is well-controlled theoretically [11] and it predicts an 
algebraic decay of the density, which follows n{t) ~ t-^/C^-i) above the upper critical 
dimension dc{k) = 2/(fc — 1) and is slowed down by fluctuations below dc{k) [11]. 
However, early simulations [23, 24] have evidenced in low dimensions the existence 
of absorbing phase transitions at non- vanishing branching rates in these systems. 
These transitions have been found to belong to two different universality classes, 
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depending on the parity of m and k §. The corresponding models are therefore caUed 
'odd' and 'even' BARW. Obviously, the effect of fluctuations cannot be neglected 
in low dimensions and need to be incorporated in a systematic way. In a seminal 
paper [12], Cardy and Tauber have derived the complete field theory for BARW 
and analysed both 'odd' and 'even' categories through perturbative RG. Their main 
results are summarised in the following. We mention here a straightforward outcome 
of their analysis, which is that for given k and m, all the reactions A: — 2, fc — 4 . . . 
and m — 2, m — 4 . . . are generated under renormalisation, the lowest m and k 
processes standing as the most relevant ones. Thus, investigating the generic processes 
[A 2A, 2A ^ 0) and {A iA, 2A ^ 0) suffices to describe the critical behaviour 
of the odd and even BARW respectively. In both cases, the general idea advocated 
in [12] is to start from the pure pair annihilation process in the vicinity of its upper 
critical dimension dc = 2 and to treat perturbatively a small branching rate a ^ e. 
The problem is then to determine if fluctuations irremediably destabilise the absorbing 
state — as suggested by mean-field — or if they allow for a non-trivial transition to 
take place. We first focus on the case of odd-BARW. 

5.1. odd-BARW: phase diagram 

For odd-BARW, the spontaneous decay to = — 1 turns out to be generated under 
renormalisation through the combination A 2A — )■ at a renormalised rate iir 
depending on the bare rates A and a. The mass A/? = cr/j — fij^ hence acquires cor- 
rections and the aim is to determine whether it can become negative for a non-zero 
bare a and drive a transition to an absorbing state. If so, the transition will naturally 
belong to the DP class since then the renormalised action is identical to (12). Cardy 
and Tauber computed A^j following two different procedures which yielded consistent 
results [12] — that are sketched on figure 2. The conclusion is the following: fluctu- 



Figure 2. Sketch of the transition lines of the BARW A ^ 2A, 2A 
computed in [12], For each dimension, the active phase Hes above the transition 
Hne, the absorbing phase below. 




ations are strong enough to induce a non-trivial phase transition only in d < 2. The 

§ Note that it has been argued in [25] that the parity conservation had no influence on similar 
reaction-diffusion systems and that the parity- wise distinction between odd and even BARW processes 
could be somewhat incidental, see also [26] . 
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transition line is given by cJc — -D[A/(2Z37re)]^/' for d < 2 and is exponentially sup- 
pressed following (Tc = £'exp(— 47r£>/A) in d = 2 [12]. The perturbation theory breaks 
down above two dimensions. However, since the annihilation rate becomes irrelevant 
in d > 2 and since moreover directed random walks are known to be reentrant only 
for transverse dimensions d <2 (that is the probability of intersection of two random 
paths vanishes above two transverse dimensions), Cardy and Tauber inferred that 
an absorbing state could no longer exist above two spatial dimensions [12]. In other 
words, their analysis suggests that the mean-field results are recovered for d > 2, that 
is the system is always in an active state above two dimensions. 

We re-examined odd-BARW using NPRG. It should be stressed that critical rates, 
like critical temperatures, are non-universal quantities. Since the NPRG formalism 
allows to keep track, through the scale integration of the flow, of the initial microscopic 
(bare) action, it enables one to determine the critical rates for which the flow leads 
to a scale-invariant effective potential. We thus integrated the flow equations dgUk, 
dgZk and dgDk at LO for different initial bare rates in dimensions 1 to 6 [2]. 

Figure 3. Phase diagrams of the BARW A ^ 2A, 2A in dimen- 
sions 1 to 6, from [2]. Lines present NPRG results, rescaled as explained in 
the text. Symbols follow from numerical simulations. For each dimension, the 
active phase lies on the left of the transition line, the absorbing phase on the right. 



d=1 d=2 d=3 d=4 d=5 d=6 




We also performed extensive numerical simulations, using an efficient algorithm 
introduced in [25] to back up our findings. For given values of X/D, we generated on 
a lattice the stochastic time evolution of particles subject to the dynamics A 2A, 
2 A ^ and we searched for the critical rates a/D which yield an algebraic decay of 
the density — in dimensions 1 to 6 [2]. 

Prior to comparing the obtained 'discrete' critical points {X/D,a/D) with the 
'analytical' ones ensuing from NPRG, the latters need to be rescaled by dimensional 
factors (C^^'', C^) to account for the continuous space limit underlying the field theory 
[2]. We emphasise that, however, matching a single point suffices to uniquely fix the 
value of the C parameter for all dimensions. Both numerical and (rescaled) analytical 
phase diagrams are displayed on figure 3 and show a remarkable agreement. Indeed, 
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let us recall that, at variance with universal quantities such as critical exponents, non- 
universal ones depend on all irrelevant operators (and microscopic details) and the 
accuracy of NPRG results was therefore unexpected. 

Let us now analyse these phase diagrams. First, the transition lines of figure 3 
are in accordance with the pcrturbativc results in their region of validity, that is in 
the vicinity of the origin. Both numerical and analytical transition lines are indeed 
quadratic in d = 1 and exponential in = 2 [2] . But away from the origin, the phase 
diagrams we obtained drastically differ from the pcrturbativc results [12] sketched on 
figure 2. Indeed, we found a transition in d = 3 and in fact in all probed dimensions up 
to = 6. The inactive phase appears to emerge only after a finite threshold (A/Z?)th 
which is in essence a non-perturbative feature. It is indeed out of reach of any pcrtur- 
bativc expansion around the origin. Furthermore, the transition lines seem to undergo 
a simple drift as the dimension grows. A closer investigation of the variation of the 
thresholds with the dimension reveals that they grow linearly with d (which has been 
checked up to d = 10) as shown on figure 4. This strongly suggests that an absorbing 

Figure 4. Evolution of the thresholds (X/D)^^ with the dimension,from [2]. 




DP transition occurs in all finite dimensions, or in other words that the mean-field 
result is recovered only in the d ^ oo limit [2] . 

We emphasise that this study unveils a remarkable instance where fluctuations 
turn out to qualitatively invalidate the mean-field and even one-loop phase diagrams. 
Indeed, fluctuations do not only bring quantitative corrections to the critical rates 
but build up genuinely non-perturbative effects that affect the very existence of the 
transition. We now come to the case of even-BARW and show that non-perturbative 
features are even more crucial there since they entirely entail the physics of the model. 

5.2. even-BARW: universality class 

The models of even-BARW have concentrated much attention [23, 12, 26] as they 
embody the first example of an absorbing phase transition not belonging to the DP 
class, but to a new one — improperly, see ^ — called 'PC for parity conserving. 
For these models, a phase transition has been observed numerically in one dimension 
and characterised by non-DP exponents. To provide an analj^tical support to these 
findings, Cardy and Tauber have attempted a KG treatment of the field theory of 
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even-BARW, which writes [12]: 

S[c^,^]^ f ${dt-DV^)(t>-\{l~^'')4>'' +a{l^$^)(t>4>. (15) 

Their analysis showed the appearance of a second upper critical dimension 
dc — 4/3 below which the branching rate a becomes irrelevant, allowing for the an- 
nihilation fixed point Fa to become stable and possibly root an absorbing state in 
d — 1. Furthermore, in a direct calculation in d = 1, Cardy and Tauber [12] managed 
to identify a combination of cr and A that admits a fixed point at one-loop order. 
However, it yielded poorly determined exponents and the extension of this calculation 
to higher orders appeared problematic [12]. 

We therefore re-analysed the even-BARW model by means of NPRG methods 
[3]. We once again exploited the generic flow equations derived in section 3, specified 
for even-BARW theory. The action (15) is Z2 symmetric. The effective potential 
Uk should hence only depend on the quadratic invariants ip'^, ^md ipipi structure 
which appears to be preserved by the flow equations dgXk's. We here consider the 
LPA, which suffices to obtain a non-trivial picture of the physics of the model. The 
effective potential is Taylor expanded around the stationary solution [tp — l^ip = 0) 
and truncated at a given power n of the fields. 



Figure 5. Variation with d of the eigenvalues of the fixed points in the 
lowest-order LPA, from [3]. Blue lines: pure annihilation fixed point F4. Red 
lines: non-perturbative fixed point FpQ whose eigenvalues are both negative for 
I < d < 1.3784 . . . and complex-conjugated at larger d (only the real part is 
plotted) . 
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Even from the lowest order n = 2 (corresponding to Uk = — Afc(l — 'i/'^)'0^ + Cfc(l — 
•f/)^ )?/;'(/)), the flow equations for the coupling constants Xk and ak exhibit, in addition 
to the Gaussian flxed point Fq = {ct = 0, A = 0} and to the pure annihilation fixed 
point Fa = {(T = 0, A^ 7^ 0}, a non-perturbative solution Fpc = {crpc 7^ 0, Xpc 7^ 0} 
[3]. The latter governs an absorbing transition in dimension d < 4/3. It becomes 
unphysical (with a negative apc) and thus plays no role for d > 4/3. Indeed, the 
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stability of Fpc and Fa can be read off from figure 5 which displays their eigenvalues. 
For d > 4/3, the annihilation fixed point is the only physical one and it is once unstable 
(in the a direction), implying that the system is always active. When d decreases, 
apc approaches zero until Fpc crosses Fa in d = 4/3, where they exchange stability. 
Below d = 4/3, Fa is hence stable, describing the absorbing phase and Fpc acquires 
an unstable direction, thus driving a phase transition [3]. The corresponding flow 
diagram is depicted for d = 1 in figure 6. 

Figure 6. Flow diagram of the lowest-order LPA in d = 1 (arrows represent 
the RG trajectories as s is decreased towards the "infra-red", macroscopic limit 
s — > — oo), from [3]. 
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The critical exponent v can be determined in the vicinity of Fpc and it converges 
as the order n of the field truncation is increased to = 2.0 ± 0.1, which is already 
in fair agreement with Monte Carlo simulations yielding v = 1.85 ± 0.1. One would 
then need to go to the next order in the derivative expansion and include non-trivial 
kinetic renormalisations analogous to Zk and Dk to get a determination of the other 
exponents. We stress that Fpc does not become Gaussian in any dimension and is 
thus genuinely non-perturbative. This explains why it cannot be reached via any 
perturbative expansion. 

6. Conclusion and prospects 

Through this review we have advocated the use of a non-perturbative field-theoretic 
method, the NPRG, to investigate non-equilibrium systems. We have highlighted the 
valuable results that this method has allowed to provide for various reaction-diffusion 
processes. Most importantly, we have unveiled that non-perturbative effects turn out 
to play a crucial role in the physics of these models. Indeed, we have first computed 
the phase diagram of odd-BARW and evidenced the existence of a non-perturbative 
threshold for the absorbing phase to emerge above two dimensions. This very threshold 
explains the failure of perturbative treatments which incorrectly lead to the conclusion 
that the transition disappears for d > 2. Furthermore, we have studied the NPRG flow 
equations of the even-BARW model and found a genuine non-perturbative fixed point 
— i.e non-Gaussian in any dimension — which is responsible for the PC absorbing 
phase transition in low dimensions and becomes unphysical above d = 4/3. It hence 
provides a theoretical back-up for the transition observed in numerical simulations. 
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reproducing in particular the numerical results for the u exponent in d = 1. 

The NPRG appears as an efficient tool to tackle non-equilibrium systems, which 
opens many prospects. First, it would allow to investigate critical properties of other 
reaction-diffusion processes. For instance, the universality class of the so-called pair 
contact process with diffusion model has been the subject of a long-lasting debate 
[27]. The influence of quenched disorder in DP models is also of great interest since 
the latter is suspected to hinder experimental realisations of this universality class 
[6]. Beyond the scope of reaction-diffusion processes, this technique has been applied 
to the study of the roughening transition in surface growth, generically modelled by 
the notorious Kardar-Parisi-Zhang equation [28]. It has allowed to obtain non-trivial 
results [29] and it would be of utmost interest to push further the investigation of 
growth phenomena. The NPRG methods could finally allow to give some insights into 
glassy dynamics. 
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